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A)  Header  information 

1 .  Name  of  principal  investigators  : 

(i)  Stdphan  Grilli,  PI  (Associate  Professor,  Dept.  Ocean  Engng.,  University  of  Rhode 
Island,  Narragansett,  RI  02882.  Tel  :  (401)  874-6636.  Fax  :  (401)  874-6837. 
Email :  grilli@mistral.oce.uri.edu) 

(ii)  Peter  Stepanishen  (Co-PI)  (Professor,  same  address) 

2.  Title  of  grant :  Modeling  of  the  acoustic  propagation  in  shallow  water  oceanic  regions 
including  effects  of  bottom  geometry  and  sub-bottom  propagation  (3/15/94  to  9/15/95) 

3.  Long  term  goal  of  the  project :  To  develop  a  numerical  tool  (model)  for  the  study  of 
acoustic  propagation  in  a  range  and  depth  dependent  shallow  water  environment,  with 
inclusion  of  bottom  effects  (scattering  and  transmission). 

4.  Category  of  research  :  propagation  (theory,  modeling  and,  computations) 

5.  Impact  on  S&T,  or  transition/integration  :  The  first  limited  numerical  model  de¬ 
veloped  in  this  project  will  make  it  possible  both  to  understand  and  to  qualitatively 
predict  acoustic  signals  measured  in  some  of  the  recent  shallow  water  experiments 
(e.g.,  Badiey  et  al,  1994).  The  new  physical  insight  gained  will  help  better  select  and 
define  the  important  features  and  the  best  approach  for  a  more  general  model,  to  be 
developed  at  a  later  stage,  that  will  address  the  long  term  goals  of  this  project. 

6.  Relationship  to  other  projects  :  This  project  has  used  results  from  recent  ONR 
supported  in-situ  experiments  of  shallow  water  acoustic  propagation. 


B)  Narrative  Documentation 
1  Background 


Emphasis  is  now  put  in  the  naval  community  towards  problems  of  acoustic  propagation  in 
shallow  water  regions.  Shallow  water  acoustic  problems  bring  new  challenges  due  to  both  the 
range  and  depth  dependent  environment  (water  stratification,  complex  bottom  topography) 
and  the  propagation  of  acoustic  waves  into  the  bottom  sediments.  For  various  reasons,  all 
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Figure  1:  Typical  configuration  of  computational  domain  for  hybrid-BEM  solution.  Domain 

1  has  arbitrary  geometry  and  is  discretized  by  BEM.  Domain  1  also  features  sources  (S)  and 
receivers  (R).  Domains  rl,  r2,  and  r3  are  radiation  domains,  in  the  water  or  in  the  bottom. 

traditional  methods  for  calculating  sound  waves  propagation  in  an  oceanic  environment  fail 
to  fully  address  these  new  challenges  (i.e.,  ray  methods,  normal  mode  methods,  parabolic 
equation  methods,  and  existing  integral  solution  methods). 

Recent  in-situ  measurements  of  low  frequency  acoustic  propagation  in  shallow  water 
over  irregular  bottom  show  significant  differences  in  the  received  signal,  depending  on  source 
frequency  and  location  (Badiey  et  al,  1994).  Most  existing  theoretical  or  numerical  methods 
for  modeling  acoustic  propagation  cannot  address  such  problems.  There  is  thus  a  need  for 
a  more  general  method  of  analysis  for  acoustic  propagation  in  a  range  and  depth  dependent 
shallow  water  environment,  with  inclusion  of  bottom  effects,  that  could  help  interpreting  and 
understanding  such  observations. 

2  Achieved  scientific  and  technological  objectives 

As  a  first  step  towards  the  long  term  goal  of  the  project,  two  simplified  two-dimensional 
models  were  developed  to  study  : 

(i)  The  effect  of  abrupt  or  gradual  changes  in  shallow  water  bottom  topography  on  the 
acoustic  propagation  in  an  homogeneous  water  column  (semi-analytical  approach). 

(ii)  The  effect  of  arbitrary  changes  in  shallow  water  bottom  topography  on  the  acoustic 
propagation  in  an  homogeneous  water  column  (numerical  Boundary  Element  approach 
(BEM);  e.g.  Brebbia,  1978). 
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Figure  2:  Typical  configuration  for  a  “step  up”  in  the  bottom. 


Approach 

A  Boundary  Integral  Equation  (BEE)  approach  was  used  for  the  representation  of  the  solution 
in  the  water  and  in  the  bottom  region  (for  porous  bottom).  Matching  was  specified  between 
the  solution  in  the  interior  region  (water)  and  the  exterior  regions  (bottom,  lateral  radiation 
regions)  along  interface  boundaries.  Eigenfunction  expansions  of  the  solution  were  used 
in  radiation  regions  (Fig.  1).  Both  analytical  (Green’s  function  and  eigenfunction  expan¬ 
sion  approaches)  and  numerical  (BEM)  solutions  were  developed  for  the  interior  problem 
equations,  with  cross-validation  between  both  approaches. 

Accomplishment  and  results 

(i)  Semi-analytical  (MATLAB)  solutions  were  developed  for  the  acoustic  propagation 
over  simple  obstacles  in  an  otherwise  flat  (rigid)  bottom.  These  were  simple  steps,  up 
and  down  (Fig.  2),  and  the  combinations  of  both  into  rectangular  obstacles,  plus  plane 
slopes  and  wedges  (Fig.  3).  Using  an  impedance  transfer  matrix,  the  method  can  also 
solve  problems  for  an  arbitrary  combination  of  steps  and  wedges  (Fig.  4). 

The  solution  of  the  problem  was  obtained  through  using  both  eigenfunction  expansion 
and/or  Green’s  function  representations  of  the  solution  in  subregions  of  the  domain 
expressed  as  an  impedance  transfer  matrix  (Jensen  et  al,  1994).  These  matrices  were 
then  combined,  using  matching  conditions  at  interfaces  between  regions,  to  represent 
the  solution  in  the  whole  domain. 

Examples  of  typical  results  obtained  using  this  approach  are  given  in  Fig.  5  and  6. 
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Figure  4:  Typical  configuration  for  a  combination  of  steps,  slopes  and  wedge  in  the  bottom. 
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Figure  5:  Example  of  amplitude  for  mode  1  for  the  propagation  over  a  step  similar  to  Fig.  2. 


Figure  6:  Example  of  amplitude  for  mode  1  for  the  propagation  in  a  wedge  similar  to  Fig.  3. 


(ii)  A  first  numerical  BEM  model  was  developed  to  address  the  same  problem  as  in  (i), 
but  for  arbitrary  (rigid)  bottom  topography  (i.e.  not  limited  to  obstacles  with  simple 
geometry). 

To  be  able,  at  a  later  stage,  to  address  problems  with  inhomogeneous  water  column, 
the  Dual  Reciprocity  BEM  approach  (DR-BEM)  was  used  to  transform  Helmholtz 
equations  into  Laplace’s  equations  (Cheng  et  al.  1994;  Grilli  et  al.  1995). 

Consistent  results  were  obtained  for  low  frequency  solutions  but  the  DR-BEM  method 
was  found  less  accurate  for  higher  frequency.  It  was  thus  decided  to  go  back  to  a  model 
which  would  still  allow  extension  to  an  inhomogeneous  water  column  but  would  be 
based  on  an  Helmholtz  equation  kernel  instead  of  a  Laplace’s  equation.  It  was  believed 
that  this  model  would  provide  better  results  for  higher  frequencies. 

(iii)  A  second  numerical  BEM  model  was  thus  developed  to  address  the  same  problem  as 
in  (i)  and  (ii),  based  on  a  direct  BIE  solution  of  the  Helmholtz  equation,  with  open 
radiation  boundaries  using  eigenfunction  expansion  representations  of  the  radiated 
fields  (Grilli  er  fl/.  1995). 

Validation  of  the  model  was  done  for  classic  analytical  solutions  with  flat  bottom 
and  the  model  was  found  to  behave  better  than  the  model  (ii),  particularly  for  higher 
frequencies. 

Many  results  have  been  obtained  following  this  approach  and  further  developments 
are  still  ongoing.  The  graduate  student  in  charge  of  this  project,  who  was  supported 
for  one  year  on  this  grant,  should  defend  his  master’s  degree  by  the  end  of  Spring  1996 
(Pedersen,  1996).  His  thesis  will  hence  represent  the  complete  exhaustive  report  on 
the  present  project  and  a  copy  of  it  will  be  sent  to  ONR  in  due  time. 
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